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Using conformal field theory, we derive several new crossing formulas at the two-dimensional 
percolation point. High-precision simulation confirms these results. Integrating them gives a unified 
derivation of Cardy's formula for the horizontal crossing probability Ylhif), Watts' formula for the 
horizontal- vertical crossing probability T\.hv{r), and Cardy's formula for the expected number of 
clusters crossing horizontally Afh{r). The main step in our approach implies the identification of the 
derivative of one primary operator with another. We present operator identities that support this 
idea and suggest the presence of additional symmetry in c = conformal field theories. 
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I. INTRODUCTION 

Percolation in two-dimensional systems remains under very active current study, despite a long history. The 2-D 
percolation point has been explored with a wide variety of methods, including conformal field theory (CFT) [1, 2], 
modular forms [3], computer simulation [2], other field-theoretic methods [4], Stochastic Lowner Evolution (SLE) 
processes [5] and other rigorous methods [6] . (We cite only a very few representative works since the literature is so 
extensive.) 

Crossing probabilities are of great interest in studies of the percolation point in two dimensions. In geometries with 
edges, these conformally invariant quantities give the probability that percolation configurations cross between some 
specified set of intervals on the boundary of the system. Perhaps the best known example is Cardy's equation for the 
horizontal crossing probability 11^ (r) [1] (which was later proven rigorously for a particular realization of percolation 
[7]). This, the probability that a percolation cluster connects the two vertical sides of a rectangle of aspect ratio 
(width/length) r, is given by 

n;,(A) = C Ai/^i^i(l/3, 2/3; 4/3; A) , (1) 

with C = 27rV3/r(i)3 = 0.56604668 . . .. The cross-ratio A is related to r by conformally mapping three consecutive 
corners of the rectangle to 1, oo and so that the fourth corner lies on the point A, with < A < 1. The interior 
of the rectangle maps to the upper half-plane. Cardy used arguments of conformal field theory; primarily that the 
(boundary) operator which changes free to fixed boundary conditions on an edge of the system is i/'i := 4>i,2 in the 
c = Kac table (the notation ipn ■= 0i,n+i simplifies the expressions for the boundary operator product expansion 
coefficients considered below [8]). 

The probability Ilhv{r) that all four sides of the rectangle arc connected by a single percolating cluster was deter- 
mined by Watts [9], using an extension of Cardy's arguments (see also the recent rigorous proof of Dubcdat [5]). This 
may be written as 

n,,,,(A) = n,,(A)-n,,-(A) , (2) 

where Ilhv denotes the probability of a horizontal crossing without a vertical crossing, 

T^mW = A 3^2(1, 1, 4/3; 5/3, 2; A) . (3) 
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To derive this result, Watts made use of a higher-order null vector in the c = = h Verma module. 

Finally, the expected number of clusters crossing horizontally, Afh{r), has also been determined by Cardy [10, 11] 
(and later via rigorous methods [12]). This calculation involves identifying percolation as the q ^ 1 limit of the 
g-state Potts model, and taking a derivative with respect to g at g = 1. Maier [13] pointed out that the result may 
be expressed as 

AA,(A) = n,(A) - in.^(A) + ^In (^^) . (4) 

The motivation for this paper is the remark by Maier [13] that the fifth-order differential equation which arises 
from the null vector used by Watts [9] to determine II^i^ has, among its additional solutions [3], both Uh and Afh- 
This mathematical observation has, to our knowledge, eluded explanation. In this paper, using a simple adaptation 
of Cardy's method, we give a unified derivation of all three formulas. In section II we calculate 11^, based on a 
physical interpretation of the V-'s 0i.4 operator. Section III extends this method to three new crossing formulas. 
Numerical verification of these results is given in section IV; by integrating them, section V reproduces the three 
known crossing quantities. Our derivation makes use of primary operators only, avoiding higher-order null vectors 
and does not require reference to the Potts models to obtain Mh- Then, in section VI, we point out that our method 
implies proportionality of 4'3 and the derivative of ipi, and explore some further consequences of this identification. 

II. CARDY'S EQUATION REVISITED 

In this section, we briefly review Cardy's derivation of the horizontal crossing probability n;i(A); and then present 
the approach used to derive it here, as an introduction to the more interesting results below. 

In [1] Cardy determined Ilh{X) via the four-point function (-0i(O)V'i(A)V'i(l)'!/'i(oo)). Here, adjacent pairs of oper- 
ators mark the intervals (0, A) and (l,oo) between which the crossing occurs. 

Our figures herein are shown as rectangles, while our formulas are given in terms of upper half-plane variables e.g. A. 
These two geometries are equivalent under a conformal mapping. One visualizes crossings in rectangles for consistency 
with common usage (e.g. "horizontal crossing" ), but takes parameters to lie on the real line for mathematical simplicity. 
Thus "the (0, A) edge of the rectangle" in fact indicates the interval on the real axis that maps into the corresponding 
side of the rectangle. Figure 1 illustrates how the four points (0, A, 1, oo) on the real axis map to the rectangle. 

Cardy's derivation (see [1] or [10] for more details) focuses on the comparison of the two possible fixed boundary 
condition assignments; either the same or different. For percolation, a fixed boundary either allows clusters to touch it 
or not; thus a rectangle with two fixed vertical edges and free horizontal edges either includes all clusters, or excludes 
horizontally crossing clusters. Therefore, by inserting a tpi (which changes the boundary condition from fixed to free) 
at each of the four corners of the rectangle, and considering the second-order differential equation implied by their 
null vector, one finds two solutions, which may be taken to be 1 and 11;, (A). Thus 

n,,(A) = (<(0)V'f"(A)V'i^(l)V'f"(oo)) - «/(0)<(A)<(l)Vf''(oo)) . (5) 

Here, the superscripts indicate the boundary condition change; / denoting free and a or b fixed boundary condi- 
tions, with a ^ b. Thus the first term includes all configurations, and is a constant, independent of A, while the 
second removes those configurations with no horizontal crossing. If we normalize our (boundary) operators so that 
{ipi{0)ipi{x)) = x~'^^\ it follows that the two sides of (5) are equal, with no multiplicative constant, and (5) becomes 

H;,(A) = 1-H;,(1-A) , (6) 

which may also be derived using duality (see [1, 10] for more details on these matters). 

Our method modifies the standard derivation as follows. Consider the probability density Wh-a that the crossing 
connects the interval (a, a -t- da) but not the interval (0, a) to the interval (1, cx)), where < a < A < 1. This is 

Hh;Q da = Hft(Q! + da) - H/i(q;) =^ Wh-a = dcJlh{a) (7) 

The configurations that will contribute to this probability are those which have a percolation cluster connecting the 
point a to the interval (1, oo) and also have, on the fixed boundary side of that percolation cluster, a dual path from 
a to (— (X),0), as illustrated in Figure 1. Here, the fact that the small interval is not connected to (0,a) ensures the 
presence of the dual path; and the differentiation in (7) removes the constant term in (5), so that a crossing cluster 
must attach to a. 




FIG. 1: Effects of (a)-configurations contributing to Tih-a- Crossing path shown in red, dual path in grey. Thick (thin) 
boundary hues represent fixed (free) edges. 



At first sight, it might seem that Ilh-.a should also depend on A. However, this quantity, which is specified in the 
half-plane, can be mapped to a rectangle with any aspect ratio r, as mentioned. When this is done, the length of the 
image of the interval (0, a) will vary according to A, which also determines r. 

Now the operator expected [6, 8] to generate a percolation cluster and dual path should have dimension h ~ I. 
This suggests that it is i/'a := 0i,4- For the moment we simply assume this, and explore its consequences. In section 
VI we give a better justification (and consider its implications). 

Note that tps, since it arises in the operator product expansion of three operators, must sit at a fixed- free 
boundary change, as shown in Figure 1. 

Therefore we have 



Uh.,o. = K (^1 (0)^3 (a)V'i (1)^1 (oo)) , (8) 

where K is a. constant. If we set 

3I/4 

it turns out that (8) will be properly normalized, as shown below. In section VI we justify (9) directly, without 
reference to percolation, by means of the operator product expansion. Note that is exactly 1/2 the constant 
appearing in (3). 
Thus 

A 



Jo 



X 



K (^.f" (0)^3"^ («)<(!)< (00)) da. (10) 



Therefore we must determine the correlation function (8). Now we may write 

■01 2^4 V'l a^3)'03(a;2)'0i a^i = 7 77 yl^ \ 1 u \] ■ 

{x:i- xi){xi- X2Y \{x:i- xi){xi~ X2) J 

Since ipi is, as mentioned, a level-two operator, the space of possible solutions for (11) is two-dimensional. However, 
in the operator product expansions of ipiipi and ipiips, the only common term is tp2 ■= 4'l,3^ so only one conformal 
block enters; i.e. the solution space is one-dimensional. To determine it, we apply the null state condition to two 
different tpi operators in (11). This gives two different second-order differential equations for F. Subtracting them so 
as to cancel the highest-order term gives 

, 2(l-2a) , , 
= F'ia) + ^^^Fia) . (12) 

This equation fixes the single conformal block as 

, (13) 
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where the superscript 2 refers to ■02 ■= 0i,3 which appears in the operator product expansion of ijji, both with itself 
and with tp^. (Our conformal blocks are normalized so that JF^" kii^) ^ x^'^~^^~^^ .) This leads to 

(V'("(0)V3^'(a)<(l)Vj^(oo)) - Ci23Cn2 (a(l - a))"'/^ , (14) 

where the usual superscripts (indicating the boundary conditions) on the boundary operator product expansion 
coefficients Cijk [14] have been suppressed, as a consequence of duality [8]. Inserting this correlation function into 
(10) then reproduces (1) 

.A 

n,,(A) = XC123C112 / a-2/3(l-a)-2/3da 

= 3XC123C112 2F1 (1/3, 2/3; 4/3; A) 

- CAi/Si^i(l/3,2/3;4/3;A), (15) 

where we have made use of (9) as weh as C123 = 2V27r/3 r[l/3]3/2 and C112 = 2,^/'^ /V[l/if/^ [8]. 

The function Tr/i;^ := dcJlh{o) is a simple example of what we call a first crossing density. The term "first" indicates 
a probability density for configurations that, when we start at the origin and move towards the point A, first contain 
a crossing cluster in the neighborhood of a. Herein, the lower case (tt) distinguishes crossing probability densities 
from crossing probabilities^ represented with upper case (H). 

With only one 03 in the correlation function we reproduce Cardy's result for 11^ (A). However by inserting an 
additional ■03 operator we can generate more complicated first (and other) crossing densities. These then give a new 
derivation of Watts' equation for Ilfiv [9], and Cardy's expression for the mean number of horizontal crossing clusters 
Nh{X) [10, 11], as weU as H,,. 



III. NEW CROSSING FORMULAS 



In order to find new results, we consider the correlation function (■0i(oo)^03(/3)03(a)^0i(O)) with < a < A, and 
I < (3. By a simple extension of the argument above, one sees that there are three configurations consistent with this 
function, illustrated in figure 2. 




FIG. 2: Configurations consistent with (?/ji(cxD)i/;3(/3)-03(a)-0i(O)). 

Let 7rJ^(a,/3) (7r^(a,/?)) denote the first crossing probability density for configurations of type A {B); first crossings 
from a to P that also make (do not make) contact with the bottom edge of the rectangle. 

Similarly, Vh{oL^(3) denotes the crossing density of configurations of type C. Now Vh{a,l3) is not a first crossing 
density; rather it includes configurations with crossings from a to (3 that are not the first crossing, but are distinct 
from previous crossings-hence the notation v in place of tt. Thus configurations with multiple crossings contribute to 
Vh{oiTp) for each pair of values a and (3 spanned by a new cluster. Integrating it therefore counts configurations with 
n horizontal crossings n — 1 times. We use this below to calculate Nh{X)- 

Now the correlation function 

(^r(x4)03(x3)^3(x2)^r(...)> - ( . y r ^.,) \ ^^^^ 

\[X4 - X2){X3 - Xi)J \{X3 - Xi)[X4 - X2) J 

It follows that 



(V'i(oo)7/.3(/?)03(a)0i(O)> = l3-^F{a/(3) 



(17) 
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Utilizing the second-order null vector for t/'i we find 



3x{l~x) ' ' 3x2(1-2;) 
Solving and selecting the appropriate conformal blocks gives 



^4 . ^ 5(l + 2x-(l-x2)2fi(l,4/3,5Ax)) 

•^13,31 - 6(1 ^^"^ 



^0 ^ _ (l + 2x+(l-x2)2Fi(l,4/3,5Al-x)) 

■^33,llU _ 3(1 ^^^^ 

•^13.11(1 -•^) = 2(1 - a05/3:,2/3 ' (22) 

with superscripts defined as in (13). The crossing symmetry relations for these conformal blocks follow using hy- 
pergeometric identities [15] for x I — x, and may be written using the operator product expansion coefficients 
C123 and C112 quoted above; we also make use of C233 = /3^^^T[l/3f^^ and C134 = a/275 [8] (note that 

K = C112/3C123, see (9)). (We have exphcitly verified that the hypergeometric identities are consistent with these 
values.) Thus 

CL-^13,3l(^) = C'n2C233^3'3,ll(l-^) (23) 
CL-^13,3l(^) = -^33,11(1 - ^) - C'n2C233^3'3,ll(l " ^) (24) 
•^33,11 (-'^ ~ ^) ~ ^123-^13,31 (^) + ^134-^13,31 (^) ■ (25) 

Given these blocks, we may use the boundary conditions to determine which configurations they correspond to. 
Fixing both intervals (0, a) and {(3, oo) in the same way determines the conformal block, so that 

{i,t{oo)i,^/mtia)ri^m = r'-^3°3,ii(i - «//?) • (26) 

With the same boundary condition on these two intervals none of the configurations in Figure 2 are excluded. Thus 
(26) is proportional to the sum of all three crossing densities. 

Multiplying by (see (9)) again results in proper normalization, as explained below. Hence 

4(a,/?) + 4(«,/3) + Ma.P) = A''/3"'-^33,ii(l " «//?) • (27) 
On the other hand, fixing the two intervals (0, a) and (/3, oo) differently leads to 

{tP(''{^)i;'/mt{a)riH0)) = Cn2C233/3-'^3'3,ii(l " «//?) • (28) 
In this case configurations of type C are excluded so that 

4Xa,f3)+7rl{a,P) = K^Cu2C233r^Tl:,,ni^ - a/ p) . (29) 
Using (24), (27) and (29) we can now find the crossing density 

Ma,P) ^ K'-p-^T°,^,,{l^a/l3)^Cn2C233ri3,ni^-a/P)) (30) 

= K^C^,,^p-'Tt,^MP) ■ (31) 

To separate 7rJ^(a, /3) and 7r^(a, f3) we fix the boundary conditions on the bottom edge (— oo, 0) to differentiate first 
crossings that touch the bottom edge (type A) and those that do not (type B or C). 
The two-point function 

{i^'/{a)^tm ={P-ar' (32) 

includes clusters connecting a and (3, but not touching the bottom edge, since it is part of a single fixed interval 
isolated by a dual path. This excludes crossings of type A, so that 

irlia, (3) + Ma, (3) = K\(3 - a)-^ . (33) 
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FIG. 3: Hull-generating walk algorithm to check for crossing densities on an 8 x 8 lattice, (left) A hull (the black curve) 
corresponding to a crossing cluster that does not touch left or bottom boundaries, and does not cross the top boundary, (right) 
Second hull walk (shown in yellow, below the first) showing that there are no lower crossing clusters. Blue segments: occupied 
bonds on the lattice. Red segments: bonds on the dual lattice, corresponding to vacant bonds on the original lattice. 



This leads to 

Collecting and simplifying these results gives the formulas 

iP + a) 2Fi(l,4/3,5/3,l-a//3)-2/3 



Vh[a,l3) 



47rv^ P^{l3-a) 
(/? + a)2fi(l,4A5/3,a//3) + 2/? 
47rV3 13'^iP-a) 
+ 2af3 ~ {f3^ - a^) 2^1(1, 4/3, 5/3, a/f3) 
47r^/3 P^{P ~ 



(34) 
(35) 

(36) 
(37) 
(38) 



These results are new, to our knowledge. They are sufficient to reproduce all three previously known crossing 
quantities, as we now proceed to demonstrate. It is interesting that only a single (2^^!) hypergeometric function 
enters. 



IV. NUMERICAL VERIFICATION 



To verify these results, we carried out simulations using hull walks on a square system, for bond percolation on the 
square lattice, where Pc = 1/2. For this system, a hull walk is a simple walk at 45° to the bonds that turns left or right 
with equal probability at each step, except when it encounters a site previously visited, in which case it always turns 
to avoid retracing its path. In this way the walk lays down the adjacent occupied and vacant bonds of a hull for the 
percolating system [16, 17]. We tested the functions Trf^{a,P), 7r^(Q;,/9), and Vh{a,P) for the half-plane transformed 
to a square system of side length 1, with a chosen to correspond to {x,y) = (0, 1/2) (the mid-point on the left-hand 
boundary), and 1 < < 00, so the corresponding point varies along the right-hand boundary. 

In the simulation, the walk was started on the left-hand side of the square at the point (0, 1/2). The requirement 
that the hull borders a first-crossing cluster starting at that point means that the hull cannot touch anywhere on the 
entire left-hand side. Similarly, if the hull crosses the top boundary, then the trial is terminated, since that event 
corresponds to the vacant bonds of the hull touching the top, preventing a horizontal crossing. Walks that touch the 
lower boundary (indicating the cluster of occupied bonds touches that boundary) were allowed to continue. Those 
walks that touch the bottom and continue to the cross to the right-hand side contribute to tt^; those crossing ones 
that don't touch the bottom were further checked for crossing clusters below them. To do this, a second hull walk was 
started from (0, 1/2), to represent the hull of the dual crossing, as in the shaded dual-lattice paths shown in Figure 5. 
If this walk intersects the bottom, there cannot be any horizontal crossing clusters below the first simulated crossing 
cluster, and the walk contributes to tt^. Otherwise, if it doesn't touch the bottom before it crosses, there must be at 
least one lower horizontal crossing, and the event contributes to Vh- 



FIG. 4: Comparison of simulation (red dots) on a 512 x 512 lattice with theory (Equations (36), (37), and (38)). 



Figure 3 shows an example of a walk on a system of 8 x 8 bonds. Here the blue circles are lattice vertices, and 
the blue edges occupied bonds on the lattice. The red edges are occupied bonds on the dual lattice, corresponding 
to vacant bonds on the original lattice. The figure on the left shows, in black, the walk corresponding to a crossing 
cluster. The system is prepared by setting one vacant bond (or a dual-lattice bond, red) immediately below the lattice 
point corresponding to (0, 1/2), so that the walk is guaranteed to enter the system and the point y = 1/2 will be at 
the boundary between occupied and vacant bonds. (Walks that exit the system at the entry point are discarded.) 
The walk then generates the remaining occupied and vacant bonds of the hull. This particular walk terminated when 
it intersected the right boimdary. Next, to check if it was a, first- crossing cluster, a second hull was initiated, starting 
on the vacant (or dual-lattice) bond in the first column. This hull is shown in yellow, in the figure on the right. 
To keep the second walk from leaving the system on the left, we added occupied bonds in the lower first column. 
This particular walk reached the bottom before reaching the right-hand side, indicating that there were no additional 
crossing clusters below the first crossing cluster. 

Because this method generates only the hull of the cluster, and simultaneously yields the type of crossing, it is very 
efficient. In several days of computer time, we were able to generate 3.3 x 10^^ hulls on a lattice of 512 x 512 bonds. 
Only 0.0016505 of the walks succeeded in making it across without hitting the top or left-hand sides. Of these, a 
fraction 0.6456 hit the bottom and contributed to tt^, while the remaining 0.3544 crossed without hitting the bottom. 

Of the latter, a fraction 0.92982 did not have additional clusters below them (contributing to tt^) and 0.07018 did 
(contributing to Uh). In all, only a fraction 0.0001707 of all initiated walks corresponded to events that contribute to 

Vh- 

The above fraction of multiple crossing events, 0.07018, is somewhat above the predicted value 0.069189, which 
is found by integrating the formulas for and ttj^ and taking the ratio of the integral of the former to the sum of 
the integrals of the former and latter. This difference can be attributed to finite-size effects, which is apparent by 
considering this quantity for lattices of side length L = 64 (0.07653), 128 (0.07294), 256 (0.07108), and 1024 (0.06978). 
The data fit very well to a straight line when plotted as a function of 1/L, with an intercept of 0.06928, quite close 
to the predicted value. 

In Figure 4 we compare the numerical results with the theory. The data are plotted vs. the position of the point 
on the right-hand side corresponding to where on the left-hand side we pick the mid-point, as mentioned above. 
The continuum coordinate was taken to be y = (y -I- l/2)/512, where the lattice coordinate F = 0, 1, . . . 511. The 
relative difference between the two curves is on the order of 2%, except near the corners of the square and where the 
functions are small, in which case the difference is somewhat larger. The overall deviation in v/i compared with the 
theory is also a finite-size effect which extrapolates nearly to zero when L — > oo. There is also a slight bias to our 
results reflecting the fact that a finite system is not perfectly symmetric with respect to the boundary conditions of 
the walk entering and leaving the system. We have found that this bias also diminishes as the system size increases. 

In conclusion, we find very good agreement between simulations and theory for these various quantities. 

V. UNIFIED DERIVATION OF CROSSING FORMULAS 

Next we integrate our formulas, to re-derive the known results for the horizontal crossing probability Ilh, the 
horizontal- vertical crossing probability Hhv, and the expected number of horizontal crossing clusters A/'/i(A). 

Now n^(A), the probability that there exists a horizontal crossing cluster that also touches the bottom edge of the 
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rectangle (such as the one illustrated in Figure 2 [A]), is given by 

/•A /"Oo 

(/3 + «) 2i^i (1, 4/3, 5/3, 1 - a/P) - 2/3 



47rV3 P^{l3-a) 



d/3 da . (39) 



(Note that there can only be one such cluster in any configuration, so n^(A) is also the expected number of this type 
of cluster.) Next let /3 ^ a/^, so that 

n^QA r ^ r (1 + 2^1(1,4/3,5/3,1-0-2 

then transform the hypergeometric function with the same identities used in deriving the crossing symmetries (23-25). 
This gives 



9a 7o (1-0^/^^^/^ 

/^r(^ + ^)-^^(/'^/f'^/^'^) + ^eda. (41) 
/o 4^V3aio (1-0 

(The coefficient of the first integral is given in terms of Cfi2 fo'^ reasons that will be clear shortly.) The identity 

3d3f,f.(M/3.5AO)= "^^'-^';;'!y^^''^'^^ (42) 

follows from the integral representation of the hypergeometric function. Using it in (41) leads to 



73 
47r 



/ 2^1 (1,4/3, 5/3, a) da. (43) 
Jo 



By (15), the first term equals n/i(A). To evaluate the second integral we use the identity 

9„ (a 3-^2(1, 1,4/3; 5/3, 2; a)) = 2-^1(1, 4/3, 5/3, a) , (44) 
which is easily derived from the scries for the hypergeometric function. The final result is 

n^,(A) = n„(A)-^A 3^^2(1, 1,4/3; 5/3, 2; A) 
47r 

= Iih{\)-Ulhv{\) , (45) 

where we have made use of (3). 

The treatment for 11^, (A), the probability of horizontal crossing when the lowest spanning cluster does not touch 
the bottom, follows analogously. Integrating 7r^(a,/3) over a and (3 as above, we arrive at the second term in (41). 
Thus 

4(A) = in;,, (A) . (46) 

Equations (45) and (46) allow us to derive Iihv{X). The configurations that contribute to n^(A) (Figure 2[B]) are 
such that the crossing from a to /3 is the first and does not touch the bottom edge. Thus the dual path from a to /3 
must itself touch the bottom edge. Therefore, by duality, 7r^(a,/3) is the probability density of a horizontal crossing 
that touches the bottom but is separated from the top by a dual cluster from a to (3. Thus n^(A) = where 
n^jjj(A) is the probability of a horizontal crossing cluster that touches the bottom, but is prevented from crossing 
vertically by a horizontal dual path. 
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FIG. 5: The five distinct configurations tiiat contribute to Xlh. Patlis in clusters are red, dual patiis grey. 
Finally, Hhvi^) is the probability of a horizontal crossing cluster that touches both the top and bottom. Hence 

n,,„(A) = niix) nL(A) = n^(A) - n^(A) 

= n^(A)-n^.B(A) . (47) 

Thus, by integrating and combining our new first crossing densities, we arrive at Watts' equation (2) for the horizontal- 
vertical crossing probability. 

Equations (45) and (46) can also be derived by a duality argument, which is a non-trivial check of our results. 
To do this, extend our notation, as shown in Figure 5. The b and t (b and t) superscripts denote configurations for 
which there is a horizontal crossing cluster which touches (does not touch) the bottom or top edge of the rectangle 
respectively. The four rightmost diagrams in Figure 5 include all the configuration types consistent with Wuv- 

Thus 



T-rbt 



bt 
hv 



jrbi 



and 

yrbt 



bt 
hv ' 



But by duality = li^^^ and = , from which (45) and (46) follow. 

Finally, we derive the expected number of horizontal crossing clusters using vuioL^ [3). Recall that this density gives 
the probability that there is a new cluster spanning from a to /3 that is not the lowest crossing cluster in the rectangle. 
Thus integrating it gives a contribution of n — 1 for each configuration with n crossing clusters. Therefore 



Vhioi, P)dl3 da 
( ^/3 



V3, 
4^^°S 



I 47r (/? 
1 

1 - A 



-4(a,/3) dp da 



1, 



- ^n„,(A) 



(48) 



giving (4). 

This concludes our derivation of the crossing formulas. As mentioned, by exploiting our new crossing results, we 
obtain all three known results without reference to the g-state Potts model or use of higher-order null vectors. Next, 
we consider our use of i/^s above from an operator point of view, and examine some of its consequences. 



VI. OPERATOR IDENTITIES 



In this section, we first consider our use of the V'a operator in sections II and III, and then present a calculation of 
the constant K used to normalize our densities (see (8), (9), and section III). 
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To begin, consider (8), which, in Ught of (5) and (7), can be interpreted as replacing dztpi{z) by Kip3{z). Now 
generally, this would not be possible, since the derivative of a primary operator is not primary itself. However the 
derivative of a primary operator of weight zero (like ^pi ) is indeed primary. 

Next, (5) gives 

= (Vi(0)i-iV'i(«)V'i(l)^i(oo)) . (49) 

Now the weight of L^itpi is 1, the same as for "03. More importantly, the null operator for -03 ;= 0i_4 is 

X>i,4 = 3Lii - 20i_2i^i + 24L2_2 + 24L_3L_i - 24L_4 . (50) 

(Here P^.s denotes the null operator for the (l>r,s Kac operator.) Further, by the Lm commutation relations (for c = 0) 
one has 

I?i,4i_i = (Lii-6L_2i-i + 6L_3)(3L?.i-2i_2), i.e. 

I?l,4l?l,l = I?3,l2?l,2 . (51) 

The right hand side is exactly the level five null operator used by Watts [9]! Since Vi ^ = 3L^i — 2i_2 is the null 
operator for ipi, so is I?i,4l?i.i as well. 

Therefore the weight of L_i^i equals that of ■03, and they both obey the same null state. Thus correlation functions 
involving them obey the same differential equations, and the solutions must overlap. Hence we posit 

L_i0i(a;) =if03(a:) . (52) 

In section VII, we discuss implications of this equation. For the moment, consider the question as to where in the 
above it actually makes a difference, i.e., if we were to differentiate a correlation function containing ipi instead of 
substituting -?sr03 for it, what would change? It is easy to see that the results of section II would be the same; however 
a crucial difference occurs for (20). Here the conformal block JF"*, which contributes to tt^, ttJ^ and Vh, would not 
appear, and our calculations would not be valid. 

Now we determine the constant K by comparing leading terms in the operator product expansions 

(L_i^i(x))Vi(o) = dMx)M0) 

= 9,(1(0) + lx^T{0) + ■■■ + ^112X^/^2(0) + . . . ) 
5 

= ^xT(O) + • • • + %ax-"2/302(O) + . . . , (53) 
6 

and 

V'3(x)0i(O) = Ci23a:-2/302(O) + • • • + Ci342:V'4(0) + . . . . (54) 

Thus 

L_iV'i(x) = T7^03(a;) = -^^3(2;) , (55) 

OL'123 ^V"" 

SO that K is indeed given by (9). Note that it appears as a ratio of boundary operator product expansion coefficients, 
rather than the derivative of the weight hi :~ ^(1,2) with the respect to the Potts parameter q, as in [10]. In fact our 
result for K also implies that 

K{i) = lih2{i)^r, (56) 



where the evaluations are at g = 1. 



VII. DISCUSSION 



In this section, we discuss a few implications of our calculations above, especially the relation (52) (see also (55)). 
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The full consequences of (52) remain to be explored. However, this relation appears to be supported by represen- 
tation theory, according to which the highest-weightspaces of a Verma module are one-dimensional [19], so that any 
two primary operators of the same weight must be proportional, as in (52). It is also interesting that the integral 
weights for the c = primary operators are exactly the Euler pentagonal numbers [20]. There are indications that 
relations similar to (52) hold for all of them. This suggest the presence of some additional symmetry for conformal 
field theory with c = 0. 

Next, consider the seventh-order null vector, which again factorizcs in two ways: 

I?3,2l?ia - I?i.52?i,2 . (57) 

Thus, arguing as above, one finds that L^iipi obeys I?3,2 as well as I?i,4 null vector conditions. 

Consider now the fusion rules of an arbitrary Kac table operator with 0i 4 and 03_2- In general one has 

[4>1a] X [(t'r.s] = [(l>r,s~3] + [0r,s-l] + [0r,s + l] + [0r,s+3] (58) 
[4>3,2] X [^r^s] = [4'r-2.s-l] + [4>r,s-l] + [0r+2,s-l] + [0r-2,s + l] 

+ [4'r,s+l] + [4>r+2.s + l] (59) 

The above then implies that only families present in both of these should be contained in the fusion rule. This 

leads to 

[L_iV'l] X [(t)r.s] = [(j)r,s-l] + [(jjr.s+l] = H'l] X [(j)r,s] ■ (60) 

Since tpi and belong to the same conformal family, they should transform among the same conformal families 

under fusions, in agreement with (60). 

Thus, our use of -03 to obtain the crossing densities augments the [ipi] conformal family. The two additional families 
present in (58) generate crossing configurations that are more complicated than those that can be generated by ipi 
operators alone. Specifically the inclusion of the '03 operator allowed us to make use of the [ipj] x [4>r,s] = [4>r,s+3] 
fusion which gives configurations of the type shown in Figure 2[C]. 

We can also use the actions of the fifth and seventh level null vectors on the identity operator to deduce properties 
of the stess tensor T. Now 

I?i,2l(z) = (3i^i -2L_2)l(z) 

= -2T{z). (61) 

Using (51) and (57) then shows that the stess tensor is annihilated by both V^^i and I?i,5. (Note that when c = 0, T 
is a primary operator.) 

Further, as argued for ip^ oc L_i'0i, only the conformal families contained in both 03_i and 0i.5 fusions should 
appear in fusions with the stress tensor, which yields 

[T] X [0,,,] = [0,,,] . (62) 

This is as expected, since the stess tensor generates conformal transformations of conformal families amongst them- 
selves. 

We hope to explore, elsewhere, the consequences of these remarks, including the "overlap" of T and 04 :~ 0i^5 in 
analogy with the result for L_i0i and "03. 

VIII. ACKNOWLEDGMENTS 

We thank A. Rocha for useful conversations. 

This work was supported in part by the National Science Foundation Grants Nos. DMR-0536927 (PK) and DMS- 
0553487 (RMZ) . 



[1] J. L. Cardy, Critical percolation infinite geometries, J. Phys. A 25 L201-206 (1992) [arXiv; hep-th/9111026]. 
[2] Peter Kleban, Jacob J. H. Simmons, and Robert M. Ziff, Anchored critical percolation clusters and 2D electrostatics, Phys. 
Rev. Lett. 97 115702 (2006) [arXiv: cond-mat/0605120. 



12 



[3] P. Kleban and Don Zagier, Crossing probabilities and modular forms, J. Stat. Phys. 113 431-454 (2003) 
[arXiv: math-ph/0209023] . 

[4] B. Duplantier, Higher conformal multifractality, J. Stat. Phys. 110 691-738 (2003) [arXiv: cond-mat/0207743] ; Conformal 

fractal geometry and boundary quantum gravity, preprint [arXiv: math-ph/0303034] . 
[5] J. Dubedat, Excursion decompositions for SLE and Watts' crossing formula, Probab. Theory Relat. Fields 134, 453-488 

(2006) (DOI: 10.1007/s00440-005-0446-3) [arXiv: math.PR/0405074] . 
[6] M. Aizenman, Scaling Limit for the Incipient Spanning Clusters, in Mathematics of Multiscale Materials; the IMA Volumes 

in Mathematics and its Applications (K. Golden, G. Grimmett, R. James, G. Milton, and P. Sen, eds.), Springer (1998) 

[arXiv: cond-mat/9611040] . 

[7] S. Smirnov, Critical percolation in the plane, C. R. Acad. Sci. Paris Ser. I Math. 333 no. 3, 239-244 (2001). 
[8] Jacob J. H. Simmons and Peter Kleban, in preparation. 

[9] G. Watts, A crossing probability for critical percolation in two dimensions, J. Phys. A; Math. Gen. 29, L363-L368 (1996) 
[arXiv: cond-mat/9603167] . 
[10] J. Cardy, Conformal invariance and percolation, preprint (2001) [arXiv: mathph/0103018] . 

[11] J. Cardy, Linking numbers for self-avoiding loops and percolation: Application to the spin quantum Hall transition, Phys. 

Rev. Lett. 84 35073510 (2000) [arXiv: cond-mat/9911457] . 
[12] S. Smirnov, unpublished. 

[13] Robert S. Maier, On crossing event formulas in critical two-dimensional percolation, J. Stat. Phys. Ill 1027-1048 (2003) 
[arXiv: math-ph/0210013] . 

[14] David C. Lewellen, Sewmg constraints for conformal field theories on surfaces with boundaries, Nuc. Phys. B 372, 654682 
(1992). 

[15] "Handbook of Mathematical Functions", M. Abramowitz and I. A. Stegun, eds.. National Bureau of Standards Applied 
Mathematics Series (1968). 

[16] R. M. Ziff, P. T. Cummings, and G Stell, Generation of percolation cluster perimeters by a random walk, J. Phys. A: Math. 
Gen. 17 3009-3017 (1984). 

[17] P. Grassberger, On the hull of two-dimensional percolation cluster perimeters, J. Phys. A: Math. Gen. 19 2675-2677 (1986). 
[18] Jacob J. H. Simmons, Peter Kleban, Kevin Dahlberg, and Robert M. Ziff, The density of critical percolation clusters 

touching the boundaries of strips and squares, preprint [arXiv: 0704.0901] . 
[19] Alvany Rocha, private communication. 

[20] Alvany Rocha-Caridi and Nolan R. Wallach, Characters of irreducible representations of the Lie algebra of vector fields on 
the circle. Invent, math. 72 57-75 (1983). 



